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We demonstrate the sensitivity of transverse-field muon spin rotation (TF-/tSR) to static charge- 
density-wave (CDW) order in the bulk of 2H-NbSe2. In the presence of CDW order the quadrupolar 
interaction of the 93 Nb nuclei with the local electric-field gradient is modified, and this in turn 
affects the magnetic dipolar coupling of the positive muon to these nuclei. For a weak magnetic field 
applied parallel to the c-axis, we observe a small enhancement of the muon depolarization rate at 
temperatures below the established CDW phase transition. Aligning the applied field perpendicular 
to the c-axis, we observe a sensivity to static CDW order in regions of the sample extending up to 
nearly 3 times the CDW transition temperature (To). The results suggest that the muon is mobile 
over the temperature range explored above the superconducting transition temperature ( T c ), and 
becomes trapped in the vicinity of defects. 

PACS numbers: 71.45.Lr, 74.70.Ad, 76.75.+i 


I. INTRODUCTION 


There is currently much interest in the interplay of 
CDW and superconducting states, driven by the discov¬ 
ery of short-range CDW correlations in the pseudogap 
phase of cuprates . 1-7 Upon lowering the temperature the 
CDW correlation volume expands, but a reversal of this 
trend below T c indicates a clear competition of the CDW 
order with superconductivity. Competing CDW order 
and superconductivity has also been demonstrated in the 
transition-metal dichalcogenide 2 ii-NbSe 2 ,— which like 
the cuprates is quasi two-dimensional. With decreasing 
temperature, 2 ff-NbSe 2 undergoes an incommensurate 
CDW transition at T 0 ~ 33 K and becomes superconduct¬ 
ing below T c ~ 7 Ki^ In contrast with cuprates, supercon¬ 
ductivity in 27J-NbSe2 does not coexist with magnetism 
or lie in close proximity to a magnetic phase, and is me¬ 
diated by a conventional electron-phonon interaction. 

The occurrence of a long-range CDW state in 2 H- 
NbSe 2 for an extensive range of temperature above 
T c differs from the situation in cuprates, where only 
short-range CDW correlations are observed above T c . 
On the other hand, 93 Nb nuclear magnetic resonance 
(NMR) measurements show pre-transitional broaden¬ 
ing of quadrupolar transitions below T ~ 75 K, which 
has been attributed to precursor CDW fluctuations 
Berthier et al. proposed that the CDW fluctuations in¬ 
ferred from the NMR measurements are induced near 
impuritieswhich bears some resemblance to a more 
recent proposal by Tacon et al. that CDW correlations 
in cuprates are pinned in nanometer-size regions by lat¬ 
tice defects^ More recently, nanodomains of static CDW 
order have been observed near defects in 2 R-NbSe 2 at 
temperatures To < T < 100 K by scanning tunneling mi¬ 
croscopy (STM)ji^ Intercalating 2 I 7 -NbSe 2 with Co or 
Mn ions reduces To, but short-range CDW order is found 
to persist for an extensive range of temperature above 
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The short-lived (with a mean lifetime of approximately 


2.2 /zs), spin S = 1 / 2 , positive muon (pi + ) is most of¬ 
ten used in condensed matter physics as a pure local 
magnetic probe, since it does not possess an electric 
quadrupole moment. However, magnetic dipolar cou¬ 
pling to quadrupolar nuclei provides an indirect sensi¬ 
tivity of the /z + to CDW correlations in the bulk. In 
2 iTNbSe 2 the spin of the implanted // + dipolar couples 
to the 93 Nb and ,7 Se nuclear magnetic moments. The 
93 Nb nuclei have spin / Nb =9/2 and possess an elec¬ 
tric quadrupole moment ( 77 Se, which has nuclear spin 
J Se = 1 / 2 , does not), which interacts with the local 
electric-field gradient (EFG). In the presence of an ap¬ 
plied static magnetic field Bo, both the /z + and 77 Se 
spins precess around B 0 , the direction of which serves 
as a suitable quantization axis. On the other hand, the 
quantization axis for the 93 Nb nuclear spins is oriented 
along the vector sum of Bo and the direction of the max¬ 
imal local EFG. 

In a TF-/uSR experiment, where the applied magnetic 
field B 0 (chosen to be along the 2 direction) is perpen¬ 
dicular to the initial muon spin polarization P(0) (chosen 
here to be along the x direction), the time evolution of the 
muon spin polarization for 2 If-NbSe 2 is well described by 

P x (t) = exp(-cr 2 t 2 ) cos(7 M B M t + </) 

= exp[-(CT^ b + a| e )t 2 ] cos^B^t + </>), ( 1 ) 


where 7 ^ is the muon gyromagnetic ratio, B M is the mean 
magnetic field sensed by the muon ensemble, (f> is a phase 
constant, and <7Nb (crg e ) is the Gaussian depolarization 
rate associated with the dipolar coupling of the // + to 
the 93 Nb ( 77 Se) nuclei. The muon depolarization rates 
are related to the width of the distribution of nuclear 
dipole fields in the 2 direction at the // + site via the 
relation cr 2 = 'y fJ ,{B 2 ), and are given by 
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-(4 Nb )j 3sin ^ cosi,j ) 2 \ . (3) 

Here 7s e (7Nb) is the gyromagnetic ratio of the 77 Se 
( 93 Nb) nuclei, 77 (rj) the distance between the /r + and 
the * th 77 Se (j th 93 Nb) nucleus, 6i ( 9j ) is the angle be¬ 
tween Bo and the straight line connecting the spin 
and the i th 77 Se (j th 93 Nb) nuclear spin, J Se = 1/2 is 
the nuclear spin of 77 Se, (lf h )j and (Ix)j are the ex¬ 
pectation values of the 93 Nb nuclear spin parallel and 
perpendicular to the z direction, and () av is the average 
over all eigenstates of the j th 93 Nb nuclear spin. Note 
that in addition to the static spin component (lf h )j par¬ 
allel to Bo, the quadrupole interaction of the 93 Nb with 
the local EFG results in a static spin component {l£ b )j 
perpendicular to Bo- The second term in Eq. m then 
represents the contribution of this perpendicular static 
spin component to the dipole field in the z direction at 
the i-i + site. Since the natural abundance of ' 7 Se is only 
~ 7.6 %, the dipolar interaction with the 93 Nb nuclei 
dominates the depolarization of the TF-/rSR spectrum. 

The ratio (Ix h )j/(lf h )j depends on the ratio of the 
quadrupolar and Zeeman interaction energies. For strong 
magnetic field the Zeeman interaction of the 93 Nb spin 
with Bo dominates and the quantization axis of the 
93 Nb nuclear spins is along the z direction. In this case 
<((4 Nb )i) 2 )av = 0 and <«4 Nb ),) 2 >av = 1(1 + l)/3, Which 
reduces Eq. (0 to an equation analogous to Eq. 0. Of 
interest here is the situation at sufficiently low fields, 
where the magnitude of the quadrupolar interaction of 
the 93 Nb nuclei with the local EFG exceeds or becomes 
comparable to the Zeeman interaction. The quadrupo¬ 
lar coupling constant for 93 Nb in 2ff-NbSe2 above Tq 
is e 2 qQ/h ~ 62 MHz^ and hence the 93 Nb quadrupo¬ 
lar frequency is vq = 3e 2 qQ / h2I(21 — 1) ~ 2.58 MHz. 
Consequently, the Zeeman and quadrupolar interaction 
energies are equivalent when the Zeeman frequency is 
uz = 2.58 MHz, which corresponds to a magnetic field of 
magnitude Bq = 27Ti/z/7Nb = 2.47 kG. The EFG at the 
93 Nb nuclear site in 2H-NbSe2 is due to the non-cubic 
ionic crystal lattice, but is modified by the presence of 
the and its screening charge, and by the development 
of CDW correlations. 

Since the p,SR technique does not require the applica¬ 
tion of a magnetic field, greater sensitivity of the fi + to 
CDW order may be achieved by removal of the Zeeman 
interaction. However, the zero-field (ZF) p.SR spectra for 
2ii-NbSe2 are rather complex. Higemoto et al. showed 
that the ZF-/iSR signals fit well to a two-component dy¬ 
namical Kubo-Toyabe function, but the physical inter¬ 
pretation of the various fit parameters are ambiguous.- 37 
Here we show via TF-p,SR measurements of 2ff-NbSe2 
that the complexity of the ZF-/zSR spectra are the re¬ 
sult of muon diffusion. Unlike the situation for zero field, 



FIG. 1: (Color online) Representative TF-pSR spectra 
recorded at a temperature near 150 K and for a magnetic 
field of (a) 25 G, and (b) 100 G applied parallel to the c axis 
of 2H-NbSe2. The solid curves through the data points are 
fits to Eq. Q. 


the TF-/rSR signals are well described by a simple single¬ 
component depolarization function, which greatly simpli¬ 
fies the interpretation of the data. 


II. EXPERIMENTAL DETAILS 

Our TF-//SR experiments were carried out on a pre¬ 
viously studied single crystal of 2 H-NbSe 2 r^ which has 
a superconducting transition temperature of T c = 7.0 K 
and an upper critical magnetic field of -B C 2 = 45 kG. The 
LAMPF and HELIOS spectrometers were used on the 
M20 surface /t + beam line at TRIUMF in Vancouver, 
Canada. The sample was mounted with the c axis par¬ 
allel to the /i + -beam linear momentum. The majority 
of the measurements were performed with the applied 
magnetic field Bo parallel to the c axis, since the Lorentz 
force exerted by a perpendicular field causes deflection of 
the muon beam. Consequently, measurements with Bo 
perpendicular to the c axis were limited to the applied 
field strength of 25 G. In all cases the initial muon spin 
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FIG. 2: (Color online) Temperature dependence of the muon 
depolarization rate a above T c for a magnetic field applied 
parallel to the c axis. The data for different magnetic field 
strengths have been vertically shifted to overlap above T = 
100 K. 


polarization P(0) was oriented perpendicular to Bo- 

The TF-/rSR spectra were fit to 

A(t) = a s P{t) + a b exp(-a b t 2 /2) cos^Bot + tp) , (4) 

where a s and a b are the initial amplitudes of the sample 
and background signals, P(t) is as defined in Eq. JTJ, and 
cr b and ip are the muon depolarization rate and phase 
constant of the background component associated with 
muons that stop elsewhere outside the sample. The fits to 
Eq. (J31) were performed assuming a b and the ratio a s /a b 
do not change with temperature. Examples of TF-/uSR 
spectra and the fits are shown in Fig. [T] 

III. RESULTS AND DISCUSSION 

Figure^] shows the temperature dependence of the de¬ 
polarization rate a above T c for different magnetic fields 
applied parallel to the c axis (Bo || c). There is some vari¬ 
ation of a with the applied field strength due to some 
uncertainty in determining the size of the background 
signal, which has a mean precession frequency close to 
that of the sample signal. Consequently, the data sets 
for different fields have been vertically shifted in Fig. [2] 
to overlap above T = 100 K, where CDW correlations are 
not expected. The ratio of the quadrupolar and Zeeman 
interaction energies of the 93 Nb nuclei for Bq = 25 G, 
100 G, 1 kG, and 4 kG are *'q/ i/ z = 99, 25, 2.5, and 0.6 
respectively. Hence the quadrupolar interaction is rele¬ 
vant to a varying degree for all magnetic fields considered 
in our experiments. 

The strong temperature dependence of a is indicative 
of thermally-activated hopping of the p + . With increas¬ 
ing temperature a first decreases, likely due to the p + 




FIG. 3: (Color online) (a) Temperature dependence of the 
muon depolarization rate a above T c for a magnetic field of 
25 G applied parallel to the c axis (This data is also shown in 
Fig-0- The solid curve through the data points is a guide to 
the eye. (b) Difference between a for Bo = 25 G and 100 G, 
and for Bo = 25 G and 1 kG. The solid curve from (a) was 
used to calculate these differences at temperatures where no 
25 G data exist. 


averaging over the distribution of fields it experiences 
during its lifetime as it moves from site to site (a situa¬ 
tion analogous to motional narrowing of the line width in 
NMR). Above T ~110 K the increase of a suggests that 
at these temperatures at least some of the p + move fast 
enough to reach sites near defects within their various 
lifetimes. Consequently, the enhanced a results from the 
muon ensemble sensing a broadened distribution of time- 
averaged magnetic fields. The plateau appearing above 
T ~ 200 K suggests that at these higher temperatures the 
p + hop rate becomes very fast, such that most of the p + 
reach the vicinity of a defect very early in their lifetime 
— where they presumably experience a broader distribu¬ 
tion of nuclear dipole fields. The muon hop rate is not 
influenced by the applied magnetic field, and hence the 
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FIG. 4: (Color online) Temperature dependence of the muon 
depolarization rate above T c for a magnetic field of Bo = 25 G 
applied parallel (cry) or perpendicular (ox) to the c axis. The 
solid red circles are the a± data multiplied by a scaling factor 
of 1.5. 

behavior of cr as a function of temperature is more or less 
the same for all fields. 

In actuality, below T ~ To there is some difference be¬ 
tween a for B 0 = 25 G and the data for higher applied 
fields. While this is noticeable in Fig. [2J it is more obvi¬ 
ous in a plot of the difference between a values. As shown 
in Fig. (31(b), the difference between cr for Bo = 25 G and 
data collected for B 0 = 100 G or 1 kG increases with de¬ 
creasing temperature below T~To. Changes in the form 
of the ZF-^iSR signal are also observed below T 0 , where 
there is no Bo to affect the 93 Nb quantization axisJi 
At Bo = 25 G, where the ratio (d^ b )/(/^ b ) ~ vq/vz is 
quite large, the quantization axis for the 93 Nb nuclear 
spin states is also solely determined by the direction of 
the maximal EFG at the nuclear site. At all tempera¬ 
tures there are contributions to the net EFG at the 93 Nb 
nuclear site from the crystal lattice and the electric field 
directed radially away from the /i + . The onset of CDW 
order below T 0 introduces an additional contribution that 
rotates the direction of the maximal EFG, hence chang¬ 
ing <TNb- Nevertheless, for this orientation of Bo the effect 
is quite weak. 

Figure @] shows the effect of changing the direction 
of the applied field to be perpendicular to the c axis 
(Bo_Lc). To maintain the transverse-field geometry, this 
also requires changing the direction of P(0) to be paral¬ 
lel (rather than perpendicular) to the c axis. Due to the 
anisotropic nature of the dipolar interaction between the 
ji + and the host nuclei, the fitted values of a are expected 
to show a dependence on the change in directions of Bq 


and P(0). In particular, the change in direction of Bq 
results in a change in the angles 9i and 0j in Eqs. m and 
m, and the change in direction of P(0) corresponds to 
a change in the direction of the n + spin. For Bo = 25 G, 
the quantization axis for the spin states rotates with 
Bo, but the quantization axis for the 93 Nb spin states 
determined solely by the net EGF does not. With the 
direction of the maximal EGF at the 93 Nb nuclear site 
fixed and continuing to define the direction of Bo as the 
z direction, the 90° rotation of B 0 changes the ratio of 
(J x Nb )/(/ 2 Nb ) in Eq. ©. 

To see the effect on a of changing the angle between 
Bo and the EFG associated with CDW order, we show 
the data for B 0 lc ( i.e. crj_) in Fig. Q] multiplied by 
a scaling factor that eliminates the anisotropy in a at 
high temperatures. Below T ~ 90 K there is a clear en¬ 
hancement of the scaled a± data above a «, indicating a 
change in the angle between the quantization axes of the 
93 Nb nuclear and n + spin states. This suggests that the 
fj, + is sensitive to static CDW order up to temperatures 
nearly 3 times the CDW phase transition temperature To 
consistent with the STM study of Ref. [lj, which as 
mentioned earlier shows nanometer-size regions of static 
CDW order in the vicinity of defects at the surface of the 
sample in the temperature range Tq<T< 3Tq. The STM 
measurements demonstrate an increase of the correlation 
length of the CDW regions as the temperature is lowered 
toward Tq. Assuming the same occurs in the bulk, with 
decreasing temperature a greater number of the mobile 
encounter regions of CDW order during their lifetime. 


IV. SUMMARY 

In the present work we have demonstrated the sen¬ 
sitivity of TF-/^iSR to static CDW order in the bulk of 
2F-NbSe2- The results support early NMR findings, and 
recent STM measurements at the surface of 2F-NbSe2 
showing the occurrence of nanodomains of static CDW 
order above To near defects. It is shown that the /x + 
is mobile at temperatures above T c , and may become 
trapped by defects during its lifetime. Above T 0 this 
mobility presumably allows a greater number of the im¬ 
planted muons to experience the short-range CDW cor¬ 
relations. 
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